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ABSTRACT 

We show that the burning of a small mass fraction X12 of 12 C in a neutron star ocean is thermally unstable 
at low accumulated masses when the ocean contains heavy ashes from the hydrogen burning rapid proton (rp) 
process. The key to early unstable ignition is the decreased thermal conductivity of a heavy element ocean. The 
instability requires accretion rates, M, in excess of one-tenth the Eddington limit when X\2 < 0. 1 . Lower M's will 
stably burn a small mass fraction of 12 C . The unstable flashes release ~ 10 42 - 10 43 ergs over hours to days and are 
likely the cause of the recently discovered large Type I bursts (so-called "superbursts") from six Galactic low-mass 
X-ray binaries. In addition to explaining the energetics, recurrence times and durations of the superbursts, these 
mixed 12 C flashes also have an M dependence of unstable burning similar to that observed. Though the instability 
is present at accretion rates « MEdd, the flashes provide less of a contrast with the accretion luminosity there, 
thus explaining why most detections are made at M » (0.1 -0.3)M£dd- Future comparisons of time dependent 
theoretical calculations with observations will provide new insights on the rp-process. 

Subject headings: accretion, accretion disks - nuclear reactions, nucleosynthesis, abundances - stars: neutron - 
X-rays: bursts 



1. INTRODUCTION 

The monitoring of galactic sources by the Wide-Field Cam- 
era on BeppoS AX and the All-Sky Monitor on the Rossi X-Ray 
Timing Explorer have detected large Type I bursts (hereafter 
referred to as "superbursts") from six accreting neutron stars 
(Cornelisse et al. 2000; Strohmayer 2000; Heise et al. 2000; 
Wijnands 2001; Kuulkers 2001). These bursts share a num- 
ber of characteristics; burst energies s=s 10 42 erg (roughly 1000 
times larger than a normal Type I burst); durations of a few 
hours and accretion rates w (0.1-0.3)MEdd (Wijnands 2001). 
The superburst recurrence time is not well known, though one 
recurrent burst has been seen 4.7 years later (Wijnands 2001). 
The light curves are similar to a Type I burst in spectral evo- 
lution, though the immediate decay timescale is about two 
hours; 1000 times longer than a typical Type I X-ray burst. 
The superbursts appear to be thermonuclear flashes from fuel 
at much larger depths than a typical Type I burst. For exam- 
ple, if the nuclear energy release is 1 MeV per accreted nucleon 
(£nuc = 10 18 erg g _1 ), the accumulated mass would be 10 24 g, 
implying a recurrence time of w 4 months for an accretion rate 
of M w O.lMEdd ~ 10 17 g s~\ much longer than a typical Type 
I burst. 

Flashes from pure 12 C layers have been discussed previously 
(Woosley & Taam 1976; Taam & Picklum 1978; Brown & 
Bildsten 1998) and might well apply to the pure helium accre- 
tor 4U 1820-30 (Strohmayer 2000). However, there are two 
difficulties with such a model for the superbursts from hydro- 
gen/helium accretors: (1) the recurrence times are very long, 
and the burst energies correspondingly large, unless a large heat 
flux from the crust of the star is present deep in the ocean; 
and, (2) our current understanding of the ashes from hydro- 
gen/helium burning does not point to accumulation of pure 12 C . 
Instead, Schatz et al. (1999, 2001) have shown that only a small 
amount of 12 C remains after all of the hydrogen and helium has 
burned via the rapid proton (rp) process (Wallace & Woosley 
1981); either from steady state burning or unstable type I X- 



ray bursting. Even though a small fraction by mass (typically 
X12 s=s 0.05-0.1), we show here that it is enough to trigger 
a thermonuclear runaway with energy comparable to the su- 
perbursts. Most important, however, is the role played by the 
heavy ashes from the rp process. These reduce the thermal con- 
ductivity enough to force a large temperature gradient in the 
ocean. This ignites the I2 C at much lower masses than pre- 
viously expected and reduces the energy and recurrence times 
to levels consistent with the superbursts. In addition, the con- 
ductive cooling times from these flashes are consistent with the 
long decay times observed. 

2. CARBON IGNITION IN A HEAVY ELEMENT OCEAN 

We presume that the ashes of hydrogen/helium burning in the 
upper atmosphere primarily consist (by mass) of a single domi- 
nant nucleus of mass Am p and charge Ze and a small amount of 
12 C with mass fractionX^. Schatz et al. (2001) recently showed 
that the onset of a closed cycle (called the SnSbTe cycle) in the 
rp-process naturally stops the increasing A trajectory once nu- 
clei of mass A w 104 are reached. Though a range of nuclei are 
made, those with masses in this region dominate the mix. For 
simplicity, we take 104 Ru as our fiducial nucleus (A = 104 and 
Z = 44) with mass fraction 1 -X12. 

The energy released when 12 C fuses (for simplicity, we pre- 
sume that the reaction makes 24 Mg) is Enuc ~ 5.6 x 10 17 erg g _1 , 
so that if X12 = 0.05, an accumulated mass of 3.6 x 10 25 g is 
needed for a 10 42 erg burst. At an accretion rate of 0.1 HiEdd 
onto a 10 km neutron star, the recurrence time would be « 10 
years (we define /«Edd = 8.8 x 10 4 gem" 2 s" 1 , the local Ed- 
dington accretion rate for solar composition). The pressure 
at the base of such a pile is P s=s 5.7 x 10 26 (gi4/2) erg cm -3 , 
where the surface gravity is g = g^lO 14 cm s" 2 . At such high 
pressures, the equation of state is that of degenerate, relativis- 

tic electrons with Fermi energy Ef = 1.9 MeV(2Z/A) 1 / 3 «- 



where p% = p/10 8 g cm 3 
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Supporting this pressure requires a 
1 .3 x 10 9 g cm" 3 and E F « 4.3 MeV. 



Though they exert little pressure, the nuclei set the ther- 
mal conductivity via electron-ion scattering. For a classical 
one component plasma with ion separation, a, defined by a 3 = 
3/4Trrij, where n,- = p/Am p , the importance of Coulomb physics 
for the ions is measured by 
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where T% = T/10 8 K. For this initial calculation we use the mi- 
crophysics of the liquid state; T < 173 (Farouki & Hamaguchi 
1993 and references therein), even when, in some instances 
(mostly when m < 0.3ra E dd) T > 173. 

Electron conduction dominates the heat transport in the liq- 
uid ocean. The heat flux is given by F = -KdT/dz, where K is 
the thermal conductivity (see Yakovlev & Urpin 1980; Itoh et 
al. 1983, Potekhin et al. 1999). ForXn < 1, we find 
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where we assume the electrons are relativistic, and take the 
Coulomb logarithm to be unity (a good approximation in the 
deep ocean; Bildsten & Cutler 1995). Integrating the heat equa- 
tion and hydrostatic balance gives the temperature contrast re- 
quired between two points (at pressure Pi and P\) to transport a 
heat flux F (Bildsten & Cutler 1995; Brown & Bildsten 1998) 
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where F E dd = m E dd(GM/R) = 8.8 x 10 gu erg cm" 2 s" 1 is the 
Eddington flux for a 1 .4 Mq, 10 km star. 

The flux in the ocean is the sum of compressional heating 
and heat flowing out from the crust, deposited there by elec- 
tron captures and pycnonuclear reactions (Haensel & Zdunik 
1990; Brown 2000). Since both contributions are proportional 
to the local accretion rate per unit area, we write F = Qm, with 
Q = 10 17 (2i7 erg g" 1 . Equation (51) then gives 
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In a careful study of the overall thermal balance of the neutron 
star, Brown (2000) found Qn sw 1 for accretion rates ;> 0. 1 mEdd 
(i.e., of the « 1 MeV per nucleon deposited in the crust, « 10% 
escapes through the surface; the rest is emitted as neutrinos 
from the inner crust and core). The compressional heating term 
contributes « 11 keV (r 8 /6) 2 (l MeV/£ F )(Z/44)(104/A) per 
nucleon, small compared to the heat flux from the crust. 

For m w 0.1 mEdd, the second term in equation (Eh is close to 
unity for a pure 12 C ocean (Z 2 /A = 3), but « 6 times larger 
for a 104 Ru ocean {Z 2 /A = 18.6) (ln(P 2 /Pi) w 8). Thus, an 
ocean of heavy ashes is hotter than that of pure 12 C at the same 
depth. This is critical, as it allows for ignition at lower col- 
umn densities, y = P/g. Figure [j] shows the temperature as a 
function of column depth for a pure 12 C ocean, and for a 12 C - 
104 Ru mixture with 12 C mass fraction X\2 = 0.1. In each case, 
we take Qn = 1, and show two accretion rates, m = 0.3 mEdd, 
and m = m E dd- The dashed line is the crystallization condition 
for 104 Ru . The lower thermal conductivity makes the 12 C - 
104 Ru oceans much hotter than the pure 12 C oceans. 

We terminate the thermal profiles in Figure n] at the depth 
at which the 12 C ignites. Brown & Bildsten (1998) showed that 



the condition for an unstable 12 C ignition in a neutron star ocean 
is set by competition with conductive cooling, so that a thermal 
runaway occurs when the local nuclear burning energy release, 
e nuc , satisfies de nac /dlnT > de coo \/dlnT , where e coo i = pKT /y 2 
is an approximation of the cooling rate (e.g., Fujimoto, Hanawa, 
& Miyaji 1981) and the nuclear reaction rates are as described 
in Brown & Bildsten (1998). For the thermal conductivity given 
by equation (Q), de coo i/dlnT = 2e CO oi, unstable ignition requires 
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where v = fl?lne nuc /<ilnr rs 26 is the temperature sensitivity 
of the fusion reaction. Ignition curves are shown in Figure p] 
for two cases. The lower dotted line is for a pure 12 C ocean 
and matches that in Brown & Bildsten (1998). The upper dot- 
ted line shows the ignition curve for a 12 C - I04 Ru mixture with 
X12 = 0.1. Because of the lower abundance, a larger column 
density is needed to ignite 12 C in a 104 Ru ocean for a given tem- 
perature. However, ignition occurs at smaller column depths in 
this case, because the I2 C - I04 Ru ocean is hotter. 

For the pure 12 C models of Figure 111 ignition occurs at 
y = 2.1 x 10 12 g cm" 2 (y = 2.4 x 10 13 g cm" 2 ) for m = m E dd 
(m = 0.3 WEdd)- Writing the burst energy Eb UKt = 4irR 2 yE nuc X\2, 
or -Ebui-st = 7 x 10 42 erg yi2X[2(R/10 km) 2 , and recurrence 
time ? recur = y/rh = 0.36 yr ynimEdd/'n), we find that pure 
12 C ignitions have recurrence times 0.8 years (29 years) and 
burst energies 1.5 x 10 43 ergs (1.7 x 10 44 ergs). These burst 
energies exceed the energy of the superbursts by at least a 
factor of ten. For the 12 C - 104 Ru models, ignition occurs at 
y = 1.0 x 10 11 g cm" 2 (y = 1.8 x 10 12 g cm" 2 ) for m = m Ed d 
(m = 0.3 /«EddX with recurrence time 13 days (2.2 years), and 
burst energy 6.9 x 10 40 ergs (1.3 x 10 42 ergs). 

3. ACCRETION RATE DEPENDENCES OF THE CARBON FLASHES 

Not only must the thermal settling solution reach the ignition 
curve, but the 12 C must also have survived to that depth. Thus, 
in order to obtain a flash, the 12 C lifetime to the fusion reaction, 
trmc ~ ^mic-XnAnuc must exceed the accumulation time, y/rh, at 
the ignition point. Combining this condition with equation (|5|) 
implies that m must exceed the critical value 
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When we evaluate this using K and the p, y relations for degen- 
erate relativistic electrons, the density cancels out, giving 
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This relation is robust and shows that large accretion rates are 
required for unstable 12 C ignition. In order to eliminate the 
temperature, we neglect the outer temperature in equation (0) 
and substitute for the temperature at the base in equation (^f. 
We then find that the explicit accretion rate dependences cancel 
and the requirement for an instability is just 
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At smaller Xi 2, the I2 C depletes and burns steadily before reach- 
ing ignition. Figure g shows the 12 C mass fraction required to 



achieve an explosion as a function of m and gn- The accre- 
tion rate dependence is mostly from the ignition pressure de- 
pendence on m. The top two panels in figure g show the burst 
energies and recurrence times for unstable models. The curves 
are forXi 2 = 0.1 and (top to bottom), g 17 = 0.5, 1.0 and 2.0. 

4. TIME DEPENDENT BEHAVIOR OF THE FLASH 

Having derived the conditions for an unstable ignition of a 
small mass fraction of 12 C , we now calculate the time depen- 
dent nature of the flash and estimate the time it takes for the 
energy to escape from the burning layer. The "one-zone" ver- 
sion of the time dependent heat equation is 
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in which we write the specific heat at constant volume Cy 
rather than constant pressure (an excellent approximation since 
Ef 3> kgT throughout the flash). In the case of predominantly 
high Z nuclei, the electron specific heat dominates that of the 
ions, so Cy = (Tr 2 ZkB/Am p )(kBT /Ef). 

Once ignited, the 12 C burning is extremely rapid and there is 
no time for cooling to occur while the I2 C burns. The fluid then 
evolves to a final temperature set by the total energy released 
J T ' CydT = J e auc dt = Xi2-Enuc- Using CV for the electrons and 
taking T{ ^> 7j gives 
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Thus, X12 > 0.01 is adequate to substantially change the tem- 
perature from that at ignition and trigger a thermal instability. If 
X\2 > 0.2-0.3, the peak temperatures would exceed 5 x 10 9 K 
and neutrino cooling would reduce the observed energetics of 
the event. 

After all the 12 C has burned and the fluid has reached 7f, the 
layer cools. Integrating equation (Jj) starting from 7} with only 
the e coo i term, we find d\nT = -dt/2t coo \, where 
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depends only on the depth in the star and the composition, and 
is temperature independent. The luminosity during the conduc- 
tive cooling phase is thus set by the temperature at the base 
of the layer; L = 47ri? 2 e C ooi} ; oc T 2 oc r f 2 exp(-f/f coo i). Thus, the 
observer should see a timescale of luminosity decay equal to 
equation (JXTj), as shown in Figure |} If we use Tf from equa- 
tion ( |1C| ) for an upper limit, we find the maximum conductive 
cooling luminosity 
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Due to our neglect of convective heat transport, this L c will be 
less than the peak luminosity of the burst. However, since con- 
vection cannot last much longer than the burning time, we ex- 
pect something similar to this simple luminosity evolution to 



occur during the burst decay, over times relevant to detection 
by flux monitoring on hourly timescales. 

Detection of these flashes should be easiest at low M's, as 
L c increases due to ignition at higher densities, whereas the ac- 
cretion luminosity, L acc , = GMM/R, is dropping. Hence, the 
contrast of the burst cooling tail is largest for low M sources, 
as shown in the bottom panel of Figure [|. Detecting one of 
these flares from a rapidly accreting Z source would require 
flux sensitivity at the 10% level on a timescale of a few hours 
and spectral sensitivity to distinguish that the flux rise is from 
extra thermal emission. Since these bursts will occur roughly 
once a month, it is an archival project to find the rare event in 
the accumulated observations. 

5. CONCLUSIONS AND FUTURE WORK 

We have shown that the ignition of a small mass fraction of 
12 C in a neutron star ocean consisting mostly of heavier ele- 
ments is thermally unstable and leads to flashes whose energy, 
recurrence times and durations are similar to the superbursts. 
The energy from these flashes takes a long time to escape the 
star, possibly explaining the persistent "offset" in the flux nearly 
a day after the superburst in 4U 1735-44 (see Figure 1 of Cor- 
nellise et al. 2000). We also predict flashes at high local ac- 
cretion rates which include the Z sources and accreting X-ray 
pulsars. In the pulsar case, some of the difficulties associated 
with explaining flares from LMC X-4 via pure 12 C burning (see 
Brown & Bildsten 1998) might well be alleviated in the mixed 
case we present here. 

Cornellise et al. (2000) reported an absence of regular type I 
bursts for seven days after the superburst. This might be due to 
thermal stabilization of the hydrogen/helium burning layers by 
the large heat flux from the cooling ashes of the 12 C burning. 
Paczynski (1983) and Bildsten (1995) have shown that lumi- 
nosities in excess of the helium burning flux (or L > L accl /100) 
will stabilize the burning, as then the thermal state of the burn- 
ing layer is independent of the local burning rates. This appears 
most likely at low accretion rates, where the conductive cooling 
is slowest, thus halting Type I bursts for ss 5f coo i ~ days after 
the superburst. 

There is still much to be done, including a time dependent 
study of the onset of the instability. Our initial integrations of 
the time-dependent code from Bildsten (1995) demonstrate a 
transition from steady-state burning at low m to unstable flashes 
at higher rates. In future work, this code will be used to more 
accurately calculate the transition accretion rate, m c with realis- 
tic mixes of 12 C and rp-process ashes. Only then can we make 
the comparisons to observations that will allow us to constrain 
the rp-process ashes. 
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FIG. 1 . — Thermal profiles and ignition curves in the deep ocean of an accreting neutron star. The ignition curves are shown as dotted lines for (bottom to top) 
Xyi = 1.0 and 0.1. The long-dashed curve shows where T = 173 for 104 Ru (the equivalent curve for 12 C lies off the right of the plot). The temperature as a function 
of column depth for a 12 C - 104 Ru mixture with Xyi = 0.1 is shown in the two upper curves, while the lower two curves are for a pure 12 C ocean. We .show two 
accretion rates, m = 0.3 riiEdd (solid), and m = riiEdd (dashed). We take Qi7 = 1. These profiles end at the depth where I2 C ignites, as defined by equation (pi). 
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FIG. 2. — The minimum mass fraction of l2 C ,X\2, needed to reach unstable ignition before depletion, as a function of accretion rate for Q17 =0.5, 1.0 and 2.0 i 
a 104 Ru ocean. 
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FIG. 3. — Burst energy, recurrence time, luminosity decay time (equation llllj) and L c (equation ||12||) as a function of m for X12 = 0.1 in i 
bottom) Q17 = 0.5, 1.0 and 2.0. The curves begin at the lowest accretion rate at which the thermal instability occurs in the 104 Ru ocean. 
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